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We have performed longitudinal magnetoresistance measurements on heavily n-doped silicon for 
donor concentrations exceeding the critical value for the metal-non-metal transition. The results 
are compared to those from a many-body theory where the donor-electrons are assumed to reside at 
the bottom of the many-valley conduction band of the host. Good qualitative agreement between 
theory and experiment is obtained. 


I. INTRODUCTION 

Magnetoresistance, the property of a material that 
its electrical resistivity changes when exposed to an 
external magnetic field, was first discovered by Lord 
Kelvin [l| . More recent discoveries are the giant magne¬ 
toresistance (GMR) ii , the colossal magnetoresistance 
(CMR)li, g, the tunnel magnetoresistance (TMR) 
the extraordinary magnetoresistance (EMR) Q, and the 
large magnetoresistance (LMR) Q. All these spectacu¬ 
lar properties were found in either magnetic systems or 
in special geometrical structures. 

Even in ordinary non-magnetic bulk materials there 
are interesting magnetoresistance effects. For all con¬ 
ducting pure single crystals it is experimentally found 
that the application of a magnetic induction B results 
in an increase of the resistivity p; the magnetoresis¬ 
tance ratio, or just magnetoresistance, defined as Ap/p = 
[p {B) — p (0)] /p (0) is positive. This general behavior of 
the crystalline state is in sharp contrast to the conduction 
properties of a number of heavily doped semiconductors 
where one observes a negative magnetoresistance. There 
are many different models in the literature trying 

to explain this anomalous behavior. They are all related 
to a model by Toyozawa [l^ where the conduction elec¬ 
trons scatter against localized spins. We refer the reader 
to the review article by Alexander and Holcomb [13 for 
the discussion of some of these models. The discussion 
is organized around a model which includes three main 
features: above the critical donor concentration, ric, the 
electrons are delocalized; above a second critical donor 


concentration, Ucf,, the Fermi level passes into the con¬ 
duction band of the host crystal; for ric < nu < ricb the 
electrons exist in a poorly understood “impurity band” 
leading to anomalous propeties. We proposed a different 
description (l^ where the donor electrons end up in the 
conduction band of the host already at the critical con¬ 
centration Uc- The anomalous properties on the metallic 
side of and close to the transition point we suggested 
were caused by many-body effects. Examples of anoma¬ 
lous behavior is that the resistivity, the heat capacity, 
and the spin susceptibility are all enhanced close to ric- 
Another example is the negative magnetoresistance that 
we treat in this work. 

Lately much research has been devoted to systems with 
positive magnetoresistance showin g a linear dependence 
on the applied magnetic field (^ . . Several mecha¬ 

nisms have been suggested to explain this behavior from 
geometrical [ 23 , classical quantum [13 113 j 

effective medium [13 perspectives. 

In this work we focus on the magnetoresistance of heav¬ 
ily doped semiconductors near and on the metallic side 
of the metal-non-metal transition. 

The material is arranged in the following way. In Sec. 
mi we present the experimental details. Sec. m is de¬ 
voted to the theoretical model and derivations. Our ex¬ 
perimental and theoretical results are compared in Sec. 
lYl Finally, Sec. |V] is a brief summary and conclusion 
section. 


II. EXPERIMENTAL DETAILS 


We use p-type, (lOO)-oriented Si wafers with resistivity 
in the range of 16-25 Hem. Phosphorous was implanted 
at room temperature. Six implantations with energies 
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Figure 1. Simulated concentration profile for nominal concen¬ 
tration 1 X 10^®cm“®. 


of 180, 120, 80, 55, 30 and 15 keV were accumulated 
in each sample with proper doses to result in a plateau 
like profile of P from the surface to the depth of 0.30/tm 
with ~ 5% deviation, according to TRIM code simula¬ 
tion [lOl- The implanted P doses were 1.4 x lO^'^cm”^ (at 
180 keV), 5.4 x 10^^cm“^ (at 120 keV), 3.6 x 10^^cm“^ 
(at 80 keV), 3.0 x 10^^cm“^ (at 55 keV), 2.2 x lO^^cm”^ 
(at 30 keV), and 1.2 x 10^^cm“^ (at 15 keV) in order to 
achieve a P concentration of 1 x 10^®cm“^. The simulated 
concentration profile is shown in Fig.[Il The doses in the 
other samples were scaled to this sample, according to the 
ratio of the desired P concentration. Samples with im¬ 
planted P concentrations of 2 x 10^®cm“^, 4 x 10^®cm“^, 
6 X 10^®cm“^, 8 X 10^®cm“^, and 1 x 10^®cm“^ were 
prepared. The damage annealing and the electrical ac¬ 
tivation of P were performed at 900 °C for 20 minutes 
in argon atmosphere in a conventional furnace. Van der 
Pauw structures [3l| were fabricated by manually applied 
indium contacts at the corners of square 6x6 mm^ sam¬ 
ples. Annealing at 300 °C on a hot plate was performed 
to improve the contacts. The implantation process as 
well as the obtained values of the electron concentrations 
are also described in detail in Refs, [s^ and (3^ . 

We performed magneto-transport measurements on 
the described structures with Van der Pauw geometries, 
exploiting conventional lock-in technique with frequen¬ 
cies 7 - 13 Hz, in the temperature range of 1.5-4.2 K and 
bias current of 10 fiA which is low enough to prevent 
heating effect and at the same time provide a well de¬ 
fined signal for our measurements. Both Hall and longi¬ 
tudinal resistance measurements were done in an Oxford 
cryostat with VTI (Variable Temperature Insert), in the 
presence of perpendicular magnetic field provided by a 
superconducting coil, capable of generating fields up to 
12 T in “^He refrigerator. 
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Figure 2. (Color online)The magnetoresistance at 4.2 K as 
a function of magnetic induction, B, for a Si:P sample with 
doping concentration 6.69 x 10^® cm~®. The red dots denote 
our experimental result and the solid curve with solid circles 
our theoretical result. The upper horizontal axis shows the 
spin-polarization parameter s. The inset shows schematically 
the energy dispersion of the spin up and spin down bands in 
presence of a magnetic field. The circles indicate which states 
are involved in the enhancement of the density of states. See 
the text for details. 

III. THEORY 

We start with our approximations and notation. Si 
is a semiconductor with v = 6 anisotropic conduction 
band valleys. For heavily n-type doped silicon, on the 
metallic side of the metal-non-metal transition (n > ric ), 
the donor electrons are up in the conduction band valleys. 
The anisotropy has some effects on the resistivity [s^ but 
we neglect this here and let the electrons be distributed 
in 1 / Fermi-spheres. The relation between the Fermi wave 
vector, fco, and the doping density, n, is given by 

ko = ■ ( 1 ) 

The Fermi energy is 

Eq = h^kl/ (2m) = k?kl/{2mdeme) (2) 

where mg is the electron mass and the density of states 

effective mass for a Fermi sphere is m^e = = 

.322. Apart from the kinetic energy there are contri¬ 
butions from the interactions between the electrons (the 
exchange and correlation energy, E^c) and from the in¬ 
teractions with the ionized-donor potentials (the band- 
structure energy, Eb). These interaction energies lead 
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Figure 3. (Color online) The same as Fig.[2]but now for 1.5 

K. 


to a deformation of the parabolic band dispersion and a 
modification of the density of states. This modification 
is important for the effects discussed in this work so we 
discuss the density of states here. 

The density of states is the number of states per energy 
and volume. The density of states from one valley is 

De = Dk/ {dE (k) /dk) = (^2n)^dE(k)/dk) /o', 

_ k'^ 

“ TT'^(dE{k)/dk) ’ 


where we have taken into account that in each valley 
there are two states for each k, one with spin up and 
one with spin down. For non-interacting electrons the 
corresponding density of states is 


D 


0 _ 
E — 


TT^ {dE° (k) /dk) 


km 

TT'^h^' 


( 4 ) 


We may express the density of states for interacting 
electrons on an analogous form by introducing a wave- 
number dependent effective mass, 


De 


km* 


The effective mass can be written as 


( 5 ) 


m* {k) = m/ [1 — /3 (fc)], (6) 


where /3 (fc) gets a contribution from each of the interac¬ 
tion energies, (3 (k) = (k) + (3b {k), where 


I3xc (k) = 
Pb (k) = 


_ rn d 5N-E^c 

TT^/c dk 5n(k) ’ 
m d 6N-Eb 

dk 5n(k) 


( 7 ) 


The quantity n (k) is the occupation number of the state 
with wave-vector k, and N is the total number of elec¬ 
trons. One effect of the interactions, that turns out to 
be very important for the present work, is that the effec¬ 
tive mass and density of states are enhanced in a region 
around the Fermi level (see Fig. 3 of Ref. 113). 

Before we continue let us introduce some dimensionless 
variables that we will use throughout. 


Q = g/2fco, 

W = nuj/4Eo, 

P = fc/2fco, 

,,^^2 


where n = 11.4 is the background dielectric constant 
of Si. In RPA (random phase approximation) the /3- 
functions at the Fermi-level are 


Pxc 


V 

Iv'TT 


1 / dQ Qg(Q 0) + TT I dQ J dW 

0 0 0 
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- 1 
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Q+1 

Q-i 


( 9 ) 

where £ (q, w) = 1 + a (q, oj) is the dielectric function 
in RPA for an electron system of v Fermi spheres em¬ 
bedded in Si with an effective dielectric constant, k. In 
the expression for j3b we have assumed that there are as 
many ionized donor potentials as there are electrons in 
the Fermi spheres and that these potentials can be repre¬ 
sented by randomly distributed pure Coulomb potentials. 

The analytical expressions for the polarizabilities 
needed in the calculations are 


a[Q,iWQ) = 


-W 


tan i^^^-btan 


W^+l-Q^ 
4Q 


w^+(i+Qy 

w^+(i-Qy^ 


( 10 ) 


and 


a(g,0) 


y 

27rQ2 



1 + Q 

[ Q 

i-Q J 


( 11 ) 


The resistivity we calculate by using the so-called gen¬ 
eralized Drude approach (36l - [^ . In the static case which 
is what we need here the results agree with the so-called 
Ziman’s formula js^- 


d a ne^Tjm* ’ 

_! 4 r 1 (1^) 

T ~ J “b:Qe2(Q^o)’ 


where p, tr, and r are the resistivity, conductivity and 
transport time, respectively. 

When a static and spatially homogeneous magnetic 
field (magnetic induction B) is applied the bands with 
spin up electrons (spin parallel to B) move up in energy 
and those with spin down electrons (spin antiparallel to 
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Figure 4. (Color online) The same as Fig.[2]but now for the 
doping concentration 6.34 x 10^® cm“®. 
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Figure 5. (Color online) The same as Fig.[3]but now for the 
doping concentration 6.34 x 10^® cm~®. 


B) move down. There is a redistribution of the electrons 
so that more electrons have spin down than have spin 
up. This has the effect that the density of states, the 
effective mass at the Fermi level, the contribution to the 
conductivity, and the transport time are no longer the 
same for the two groups of electron. Let us introduce the 
spin-polarization parameter, s, that varies from zero in 
absence of B to 1 at full polarization (all electrons have 
spin down), 


s = 


n X —n t 


n 


(13) 


The density and Fermi wave-number of spin up and 
down electrons are 


n t= 
n ],= 

ko t= fco/a, 

ko ),= ko/b, 

where 

a = (1 - 
6 = (1 -I- 

The resistivity is now 

~ nte^Tt(l-/3t) I n4-e^T4.(l-3J.) 

_ _ m/e^ _ 

“ ratTt(l-/9t)+"iTi(l-/34) ■ 


(14) 


(15) 


( 16 ) 


The polarizability can be divided into a contribution 
from each group of electrons 

a (Q, iWQ) = a t (Q, iWQ) + a ; (Q, iWQ ), 
a t (Q, flLQ) = (aQ,m^VFQ) , (17) 

a ], {Q,iWQ) = -^a {bQ,ib'^WQ) , 

and the new transport times become 


1 4 1 r 1 

rt “ 3 771^ (1-s) J Q[l + {a/2)a{Qfi) + {b/2)a{{b/a)Q.,0)\-‘ ’ 

1 4 ve^m 1 r 4/n 1 

rj. - 3¥I^IT+3) J Q[l+(b/2)a(Q,0) + (a/2)a((a/b)Q,0)Y ' 

(18) 

We also need the modified /3-functions. For spin up 
electrons we have 


’t= 


ya 

L'TT 


1 I Q[l+{a/2)a{Q,0) + {b/2)a{{b/a)Q,0)] 
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2(Q+1) 
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Q^(Q^-l)[l+(a/2)a(Q,0) + (b/2)a((b/a)Q.,0)f 

X 


1 Q In 
^ 2Q 


Q+1 

Q-i 


(19) 

To get P P one just interchanges a and b in Eq. (IT^ . 
Now we have all formalism needed for the calculation 


of the magnetoresistance, Ap/p = [p (s) — p (0)]/p (0), 
as a function of spin polarization s. The experiments are 
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Figure 6. (Color online) The same as Fig. [3] but now for a 
sample with two donors, P and Bi, with the total doping 
concentration 8.4 x 10^® cm“®. 


not given as a function of s but as a function of B. If the 
fields are small enough we can assume a linear relation 
between B and s. It can be written as 


B[T] = 


2.099879 X [cm-^] 

rride ix/Xo) 


( 20 ) 


IV. EXPERIMENTAL AND THEORETICAL 
RESULTS 

Our theoretical and experimental results are compared 
in Figs.[2] - [6l We have adjusted the spin-susceptibility 
enhancement-factor (x/xo) appearing in Eq. (|2n|l to get 
a reasonable fit between the theoretical and experimen¬ 
tal curves. The adjustment only affects the theoret¬ 
ical curves in the horizontal direction. The enhance¬ 
ment factor is the only fitting parameter we have intro¬ 
duced. The enhancement of the density of states at the 
Fermi level increases when the density comes closer to tic 
{ric = 3.5 —4.4 X 10^®cm“^ [s^) from the metallic side. It 
is furthermore well known that the spin susceptibility y 
is more and more enhanced the closer to ric one gets and 
that the enhancement is reduced when the temperature 
goes up (diil - li^ . Our extracted (x/xo) comply with this 
behavior. In the figures the upper horizontal axes show 
the spin polarization parameter, s, defined in Eq. (1131) . 
The relation between s (upper horizontal axis) and B 
(lower horizontal axis) is given by Eq. 

Now, what causes the negative magnetoresistance? As 
we mentioned above the density of states is enhanced 


at the Fermi-level in absence of a magnetic field. This 
leads to an enhancement of the resistivity. In absence 
of magnetic field the spin up and spin down bands are 
degenerate and the Fermi wave-numbers are the same for 
both spin types. When the magnetic field is introduced 
the spin down bands move down in energy and the spin 
up bands move up. There is a redistribution of electrons 
from the spin up bands to the spin down bands so that 
the Fermi-level is the same in all bands. The Fermi wave- 
numbers are now different in the two band types. See the 
inset of Fig. [21 When the magnetic field is introduced the 
density of states of both electron types, i.e. spin up and 
spin down electrons, are enhanced for states with wave- 
number ko't and ko This means that the peak in the 
density of states at the Fermi-level is for each spin type 
split up into two. The states involved in the enhancement 
of the density of states are indicated by circles in the inset 
of Fig. 121 For spin up electrons one peak remains at the 
Fermi-level while the other moves up into the unoccupied 
part of the bands. For spin down electrons one peak re¬ 
mains at the Fermi-level and one moves further down in 
the occupied part of the bands. For both spin types the 
enhancement at the Fermi-level is hence reduced. It is 
only the enhancement at the Fermi level that effects the 
resistivity. This causes the initial negative magnetoresis¬ 
tance. The enhancement of the density of states at the 
Fermi-level for both spin types as function of magnetic 
induction is shown in Fig. [71 for the sample with doping 
concentration 6.69 x 10^® cm“® at 4.2 K. 

There is another effect that acts in the same direc¬ 
tion. There are Friedel oscillations in the screening- 
charge density centered around each impurity potential 
with Fourier component q = 2fco or Q = 1. This leads 
to an enhanced scattering rate in the back scattering di¬ 
rection across the Fermi spheres and an enhancement of 
the resistivity. At zero magnetic held the Friedel oscil¬ 
lations have the same periodicity for spin up and spin 
down electrons; both type of electrons scatter equally 
strongly against both Friedel oscillations. When the mag¬ 
netic held is turned on the Friedel oscillations will split 
up into two; one with Fourier component q = 2ko f; one 
with Fourier component q = 2fco This means that 
the back-scattering rate for an electron of a certain spin 
against the Friedel oscillations of the opposite spin is re¬ 
duced. The enhancement of the resistivity is thus reduced 
leading to a negative magnetoresistance. However, this 
effect is expected to give a much smaller contribution to 
the negative magnetoresistance than the density of states 
effect since the electron can scatter with a wave number 
ranging from zero up to two times the Fermi wave num¬ 
ber. 

All our calculations are for zero temperature. What 
happens for non-zero temperatures? If we study classical 
experiments we hnd that the negative magnetoresis¬ 
tance effect is gradually reduced when the temperature 
is enhanced. This is consistent with our theory. The 
peak at the Fermi-level of the density of states is ex¬ 
pected to be broadened. Besides, at zero temperature 
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Figure 7. (Color online) The enhancement of the density of 
states at the Fermi-level for both spin types as function of 
magnetic induction, B. The results are valid for the doping 
concentration 6.69 x cm“® at 4.2 K. 


Figure 8. (Color online) The modulus of the maximum neg¬ 
ative magnetoresistance from our theoretical calculations. It 
shows linear relation on a log-log plot which means a power 
law dependence on the doping concentration . 


only states at the Fermi-level takes part in the conduc¬ 
tivity. When the temperature goes up also states away 
from the Fermi-level where the enhancement in the den¬ 
sity of states is weaker take part. Both these effects are 
expected to gradually remove the negative magnetoresis¬ 
tance. The temperature effects are expected to be more 
and more important the lower the density. This is also 
what is found experimentally . 

In Fig. [5] we see that the maximum negative magne¬ 
toresistance increases linearly on a log-log plot when the 
density is reduced. For finite temperature the maximum 
is expected to start decreasing at a density that depends 
on the temperature. The higher the temperature the ear¬ 
lier the decrease is expected to set in. In Fig. 8 of Ref. [i| 
this decrease is observed. 


V. SUMMARY AND CONCLUSIONS 

We have performed magnetoresistance measurements 
of heavily phosphorous doped silicon and compared the 
results to theory. The resistance was calculated using 
the so-called generalized Drude approach taking many- 
body effects into account. We propose that the origin of 
the negative magnetoresistance is a combination of two 
effects. The many-body effects lead to an enhancement 
of the density of states at the Fermi level which in turn 
results in an enhancement of the resistivity. Friedel oscil¬ 
lations in the screening charge density cause an enhanced 
back-scattering rate across the Fermi volumes leading to 
an additional enhancement of the resistivity. 

When the magnetic field is turned on the enhancement 


of the density of states for each spin type is split up into 
two peaks, one at the Fermi-level and one that moves 
away from the Fermi-level with enhanced magnetic field. 
This reduces the resistivity. Also, the back-scattering 
rate against the Friedel oscillations will be reduced in 
the presence of a magnetic field. Both these effects act 
towards a negative magnetoresistance. The first effect 
is expected to be dominating. We would have liked to 
have samples that cover a broader range of doping con¬ 
centrations from closer to Uc to a much higher value. 
Unfortunately, we have had contact problems preventing 
high quality measurements. We hope to come back with 
complementary results in the near future. 

From our work we can deduce that the metal-non- 
metal transition in heavily n-doped Si is driven by 
electron-electron and electron doping-ion interactions 
and not by spin-flip scattering or spin-orbit scattering. 
This is in line with the findings in Ref. for p-doped 
Si (Si:B) where these conclusions could be drawn from 
the universal scaling of the magnetoconductance of that 
system. 

We have not been able to hnd any recent work on the 
magnetoresistance of heavily doped Si near the metal- 
non-metal transition. We hope that the present work 
will inspire the readers to perform new such measure¬ 
ments. Recently [i^ one has managed to produce sam¬ 
ples of sulfur-doped silicon at non-equilibrium concentra¬ 
tions in a range covering the concentration of the metal- 
non-metal transition. This was achieved by using ion 
implantation followed by pulsed-laser melting and rapid 
resolidification. Sulfur is a deep level impurity in crys¬ 
talline silicon and tic is much higher than for phospho- 
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rus doping. One estimated Uc to be between 1.8 and 
4.3 X 10^°cm“^, i.e., two orders of magnitude higher than 
for phosphorous doping. It would be very interesting to 
see magnetoresistance measurements being performed on 
these samples. 
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